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The entanglement entropy of a gapless fermion subsystem coupled to a gapless bulk by a "weak 
link" is considered. It is demonstrated numerically that each independent weak link contributes an 
entropy proportional to lnL, where L is linear dimension of the subsystem. 



PACS numbers: 71.10.-w, 03.67.-a 



INTRODUCTION 



Striking connections have recently emerged between 
quantum entanglement and several distinct fields of 
physics: black hole quantum mechanics, quantum phase 
transitions in condensed matter, topologically ordered 
phases in quantum field theories (QFT) and conformal 
field theory (CFT). The central quantity common to 
these studies is the entanglement entropy, computed for 
a finite subregion of a QFT or many body system. If, 
in a QFT in d spatial dimensions, a distinguished re- 
gion A of volume L d is formed, it follows that the de- 
grees of freedom which reside exclusively in the region 
A will appear to be in a mixed state. The degree of 
mixing may be characterized by the von Neumann en- 
tropy, S = — trplnp, where the reduced density matrix 
p = tr^ ^4 1 0) (0 1 has been formed by tracing over the de- 
grees of freedom of the ground state, |0), exterior to the 
region A. 

Entanglement entropy was first investigated in the con- 
text of black hole quantum mechanics and Hawking- 
Bekenstein entropy where it was found that entan- 
glement entropy is not a conventional extensive quantity 
but, rather, is proportional to the area of the bound- 
ing surface, S oc [i]. This highly suggestive result 
(known as the area law) is believed to bear some rela- 
tion to holographic principle proposals ft has also 
recently been discovered that the subleading term in the 
entanglement entropy carries information about topolog- 
ically ordered phases occurringin topologic QFTs [4j . As 
pointed out in references [!, 0,0, 0] , the divergence of en- 
tanglement entropy at quantum critical points may also 
be exploited to identify quantum critical phenomena and 
quantum phase transitions. In connection with entan- 
glement entropy, 1 + 1-dimensional CFTs have received 
special attention in the past several years. CFTs — which 
describe critical spin chains, Luttinger liquids and other 
massless theories — have pointlike bounding surfaces; re- 
markably, the entanglement entropy was shown to de- 
pend universally upon the central charge of the theory 
and to diverge logarithmically with the length of the 
subsystem [J Specifically, the entropy is given by 
S — | InL/e where c is the central charge. 

So far, studies have concentrated on entanglement en- 



tropies computed in a homogeneous system; specifically, 
the bipartite entropy of a d-dimensional distinguished re- 
gion separated from the d-dimensional bulk by a (d— 1)- 
dimensional boundary. For gapless fermions it has been 
proven that the area law is anomalous and, in contrast 
to the area law above, 



(1) 



where e is a spatial cut-off [11|, |12|, [13, ll4j • One might 
consider a generalization in which a (d — fc)-dimensional 
boundary separates the distinguished region from the 
bulk, where the codimension of the boundary, k, is some- 
thing other than k — 1. To maintain the appropriate ex- 
tensivity of the entanglement entropy, one might expect 
S ~ L d ~ k \nL, although such a property has not been 
established. In this paper we study the simplest inter- 
esting subcase: a zero dimensional boundary — or "weak 
link" — connecting two d-dimensional regions (see fig. 1). 
For d-dimensional gapless fermions, this geometry may 
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1: Two d- dimensional regions connected by a "weak 



be treated by bosonization and mapped to the problem 
of entanglement in a 1-dimensional subsystem where the 
entropy is known. The treatment of the weak link is sim- 
ilar to the x-ray edge or Kondo problem, where the impu- 
rity interacts with the s-wave sector of the bulk fermion 
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system [h|. The weak link leads to the entanglement of 
effectively two (length L) 1-d systems and the entropy is 
thus proportional to InL. Each independent weak link 
contributes an entropy proportional to lni, consistent 
with the bulk area law (eq. [I} if one considers a large 
number of weak links proportional to L^ 1 . In a sense, 
the area law (eq. [I} may be anticipated by the fusion of 
two d-dimensional regions in this fashion. 



MODEL HAMILTONIAN AND BOSONIZATION 

In this manuscript, we study the entanglement entropy 
of a 2-d or 3-d noninteracting fermion system (A) con- 
nected by a weak link to a second identical system (B). 
The model hamiltonian is: 



H = -t ( C L^>H 

(i,j)a=A,B 

- 2/(cq,a c o,s + cI b c ,a) 



«) 



(2) 



where i and j are two dimensional site indices and A 
and B denote the two identical systems coupled through 
a hopping amplitude, y. For the numerical part of this 
study, each subsystem is taken to be a periodic square or 
cubic lattice of linear size L. Figure 2 shows a sample cal- 
culation of the energy eigenvalues of the hamiltonian . 
Unlike the comparable coupled 1-d systems for which en- 
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of the hamiltonian ([2}. This model, as well as its d- 
dimensional counterparts, may be bosonized in a stan- 
dard way, analogousl y to the kondo impurity model and 
x-ray edge problem 1 101 ] . Referring to figure 1, the s- 
wave in/out fermion modes (denoted R and L) may be 
combined into a single spinor 



ip%(x) x > 
-ipt i-x) x < 



where a refers to subsystem A or B. Since the fermion 
spinor now obeys periodic boundary conditions, \I> a (L) = 
^ a {~ L), the model may be bosonized in the standard 
way, 



(3) 



where (j> a (x) is a right-moving boson field on the whole 
interval [—L,L]. The hamiltonian written in chiral form 
is now 

H = j L [{d x 4> A {x)f + {d x <l> B {x)f]dx (4) 



- ycosV47r(0 A (O)-^(O)) 

This model is equivalent to the model studied analyti- 
cally at weak coupling in reference I15II and numerically 
for arbitrary coupling in reference [1 €31 ] - Alth oug h these 
studies reach similar conclusions, reference [16| defini- 
tively shows that when an impurity is introduced into a 

1- d system of gapless fermions, the entanglement entropy 
of a subsystem of length L, in which the impurity lies on 
the boundary, is proportional to In L with a nonuniversal 
prefactor depending upon the impurity coupling. From 
the above reduction to a 1-d problem, it is plausible that 
each sufficiently separated weak link connecting the two 

2- d subsystems would contribute an entropy proportional 
to In L. We expect entanglement to behave as 



S = a{y)N\nL + b 



(5) 



where the prefactor a(y) depends upon the weak link 
coupling constant, y, N is the number of independent 
weak links and b is a constant. 



FIG. 2: Eigenvalues of hamiltonian ([2]) for an L 
lattice. 
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tanglement entropy has been studied in refs. [II, [II [ill, 
the coupled 2-d systems have extensive (in L) degenera- 
cies. Moreover, the weak link only partially lifts these 
degeneracies in each shell to form two localized states. 
These features will give rise to a special dependence of 
the entanglement entro py upon the filling fraction of the 
lattice. References ll|, [l2|, ]l3|, ]l4j have closely studied 
the role of critical points in the fermi surface in connec- 
tion with the anomalous fermion area law. 

We anticipate the central numerical results of the next 
section showing S oc In L by appealing to bosonization 



COMPUTATION OF ENTROPY 

We compute the entanglement entropy of subsystem A 
following the method introduced by Peschel in ref. [l8j |. 
The entropy may be computed from the eigenvalues of 
the ground state free fermion correlation matrix 



C x ,y — (c y c x ) 



(0) 



where x and y are lattice points exclusively within sub- 
system A. Denoting by the eigenvalues of C x<y , the 
expression for the entanglement entropy is: 
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Figure 3 shows an example of the correlation function 
eigenvalues which contribute appreciably to the entropy; 
almost all eigenvalues are essentially zero or one. As first 
pointed out by Peschel [is| . the set of which may 
be thought of as an "effective" single particle distribu- 
tion function for subsystem A, have the characteristics of 
a thermal fcrmi distribution at a fictitious temperature 
(depending upon L and y in the present case.) Figure 3 
shows a fit to the fcrmi function. 
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FIG. 3: Eigenvalues of correlation function matrix (eq. [6]) for 
an L = 25 lattice with y = 1. The eigenvalues are seen to fall 
upon a finite temperature fermi distribution (solid line). 

Figure 4 shows the entanglement entropy of a single 2-d 
subsystem with linear size L coupled to an identical sub- 
system. The entire coupled system has a fermion filling 
fraction set to v = 1/2 (that is, the number of fermions 
is L 2 ); for even L, this sets the fermi level in the middle 
of a degenerate shell (see fig. 2). Since the degeneracy at 
the fermi surface is proportional to L, the entanglement 
of these states generically leads to an entropy propor- 
tional to L, as seen for the even L lattices in figure 4. 
For odd L lattices at v — 1/2, this extensive degener- 

25 -i 1 

-I 1 1 1 1 r—^ 1 1 1 1 ,— 

10 12 14 16 18 20 22 24 26 28 30 32 

L 

FIG. 4: Dependence of entanglement entropy on 2-d system 
linear size, L, for weak link hopping amplitude of y = 1. 

acy is avoided and the system better approximates the 
continuum limit leading to the bosonized action, eq. Q. 



Returning to figure 3 showing the spectrum of eigenval- 
ues of the correlation function, typical for an odd lattice 
where the extensive degeneracy is avoided, it is seen that 
the quasi-fermi level is at (L 2 + f)/2 = 313 and only 
a few eigenvalues contribute significantly to the entropy. 
Now, expanding the vertical scale of figure 4 and plotting 
only odd L data versus InL, figure 5 clearly shows en- 
tanglement entropy proportional to InL, confirming eq. 
(for N = 1) even for relatively small lattices. 



CO 1.1 




2.3 2.5 2.7 2.9 3.1 3.3 3.5 3.7 3.9 4.1 

Ln L (odd L) 



FIG. 5: Dependence of entanglement entropy on 2-d system 
linear size for odd L and weak link hopping amplitude of 
y = 1. Linear system size in this figure ranges from L — 11 
to L = 57. 

We have also investigated the entanglement entropy in 
a 3-dimensional geometry. Two 3-d subsystems (obey- 
ing periodic boundary conditions) were coupled together 
by a single weak link as described by hamiltonian @. 
The degeneracies in this system are more complicated 
than in 2-d and as a result, the degenerate shell close to 
zero energy is rarely completely filled for half-filled lat- 
tices (v — 1/2) with different linear sizes, L. Instead 
of fixing the value v = 1/2, the closest filling fraction 
to v = 1/2 was used that filled the degenerate zero en- 
ergy shell. Theses filling fractions fell within the range of 
v = 0.497—0.565 for the lattice sizes considered. Figure 6 
shows the results for 3-d entanglement entropy. Entropy, 
S, is monotonically increasing with linear system size, L, 
and exhibits a periodic variation resulting from the vari- 
ation in v required to achieve the filled shell condition 
(see inset). Although it is impossible with our present 
numerics to examine the entropy over an increase in L 
by a factor of more than approximately 3, the behavior 
of entropy looks approximately proportional to ln L for 
the lattice sizes considered. 

Lastly, we turn to the numerical estimate of the co- 
efficient a(y), studied in 2-dimcnsional geometries. The 
dependence of the prefactor in eq. (O, a(y), upon the 
weak link coupling constant, y, seems to be approxi- 
mately quadratic for small y, and is shown in figure 7. 
Similar behavior was found for 1-D systems in reference 
0. 
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FIG. 6: Dependence of entanglement entropy on 3-d system 
linear size for weak link hopping amplitude of y — 1. Filling 
factor, v ~ 0.5, is adjusted for each linear system size, L, 
to completely fill the degenerate shell closest to zero energy. 
Inset: entropies plotted without the logarithmic scale show 
periodic variation in entropy corresponding to the filled shell 
condition in 3-d. Linear system size in this figure ranges from 
L = 6 to L = 19. 




FIG. 7: Dependence of prefactor a(y) in eq. ((5J on coupling 
constant y for system sizes L — 15, 19 . Small y behavior is 
approximately quadratic. 



MULTIPLE WEAK LINKS 

We now consider the entanglement entropy of large 
lattices coupled by more than one weak link. The hamil- 
tonian is the same as eq. ^ , but modified by additional 
weak link interaction terms at different sites. Since each 
weak link lifts the degeneracy of two more states from the 
fermi surface, the number of fermions, N, must be set to 
be N = L 2 — W+ 1, where W is the number of weak links 
(and L is odd) to maintain the "closed shell" condition. 
The two-fold lifted degeneracy found for a single weak 
link may be seen in figure 2. 

Figure 8 shows that the entropy continues to be pro- 
portional to InL for two weak links, consistent with eq. 
([5]). The entropies shown for two weak links in figure 
8 were computed with the weak links separated maxi- 
mally at (x, y) locations on the lattice corresponding to 
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FIG. 8: Dependence of entanglement entropy for one and two 
weak links on 2-d system linear size for odd L. Consistent 
with zero dimensional area law, eq. 10, slope of graph for 
two weak links is double that of single weak link. Weak link 
hopping amplitude is y — 1. 



(L/4, L/4) and (3L/4, 3L/4). However, the same com- 
putation with two weak links on nearest neighbor sites 
generates results that are virtually indistinguishable from 
those of figure 8. A linear least squares fit to the S versus 
InL data in figure 8, yields a slope of 0.338 for one weak 
link and double that value, 0.676, for two weak links. We 
note that the intercept, b, in eq. ([5]) was determined to be 
0.168 and 0.316 for one and two weak links, respectively, 
possibly following a linear relationship b oc N. Such a 
behavior would not be surprising — reflecting a subdom- 
inant O(L) contribution to the bulk area law, eq. (JTJ . 
Thus the entropy accorded each weak link appears to be 
additive, largely independent of the spacing of the two 
weak links. 
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FIG. 9: Dependence of entanglement entropy upon the sepa- 
ration between two weak links (with y = 1) coupling subsys- 
tems A and B. The zero separation data point corresponds 
to the entropy of a single weak link (with y = 1) for reference. 

The separation of two weak links is explored further 
in figure 9, where the entropy is computed for two weak 
links on an L = 15 lattice, with a separation that ranges 
from to 6 lattice spacings. Other than a slight (not 
understood) odd/even modulation, the entropy is inde- 
pendent of the spacing, and its value, S ~ 2.16 — 2.19 is 
approximately double the entropy of a single weak link 
with identical coupling, S ~ 1.09. 
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To explore the additivity of the entanglement entropy 
for multiple weak links, we computed entropy for up to 
four weak links on an L = 25 lattice, maintaining the 
maximum possible separation between them. As shown 
in figure 10, each additional identical weak link con- 
tributes an entropy of approximately S ~ 1.26, until 
some saturation appears for four weak links. In terms 
of the Fermi wavelength, Xp, saturation begins to ap- 
pear at a linear separation of approximately 12 lattice 
sites, corresponding to 6Af for a half filled lattice. It is 
interesting to note that although saturation effects are 
seen at a spacing much larger than A^, two weak links 
separated by only a single lattice spacing appear to be 
independent in that they contribute entropies additively 
as shown in figure 9. 



was, in part, completed; this research was supported in 
part by the National Science Foundation under Grant 
No. PHY05-51164. 
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FIG. 11: Two 2-dimensional regions, A and B, connected 
by several weak links (shown as short radial line segments). 
Taking the boundary length to be L and choosing a number 
of weak links, N = is equivalent in the zero dimensional 
area law (eq. [5]l to choosing a cut-off e in the bulk area law 
(eq. p. 



FIG. 10: Dependence of entanglement entropy upon the num- 
ber of weak links. The weak links were maximally separated 
in an L — 25 lattice. 



CONCLUSION 

We have demonstrated that a weak link coupling two 
gapless fermion subsystems results in an entanglement 
entropy between the subsystems proportional to In L 
where L is the linear dimension of the subsystem (eq. 
03). This result is consistent with the bulk area law for 
fermions (eq. EQ) in that the entropy in a geometry con- 
taining multiple, sufficiently separated, weak links is pro- 
portional to both In L and the number of weak links. The 
consistency may also be seen by noting that the bulk area 
law depends upon an arbitrary cut-off, e, depicted in fig- 
ure 11 as a spatial cut-off taken tranverse to the bound- 
ary. Choosing a cut-off in the bulk area law (eq. [TJ, as 
shown in figure 11, is equivalent to choosing a number of 
weak links, N, in the zero dimensional area law (eq. 03) 
connecting the two regions satisfying, 

This research was supported by an award from Re- 
search Corporation CC6535. GL also wishes to thank the 
Kavli Institute for Theoretical Physics where this work 
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